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Abstract 

We have generalized the properties with the tensor product, of one 4x4 
matrix which is a permutation matrix, and we call a tensor commutation 
matrix 2 £g) 2. Tensor commutation matrices n <g> p can be construct with 
or without calculus. A formula allows us to construct a tensor permutation 
matrix, which is a generalization of tensor commutation matrix, has been 
established. The expression of an element of a tensor commutation matrix 
n <S> p, with n = p , has been generalized in the case of any element of a 
tensor permutation matrix. The tensor commutation matrix 3 (g) 3 has been 
expressed by using the Gell-Mann matrices. 



INTRODUCTION 

When we had worked on RAOELINA ANDRIAMBOLOLONA idea on 
the using tensor product in Dirac equation pQ, [2] we had met the unitary 
matrix 
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which has the following properties : for any unicolumns and two rows matrices 
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G M 2xl (K), [0] 



e M 2xl (K) 
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and [1. for any two 2x2 square matrices, [A], [B] G.M2x2(lK) 

[U 2m ]-{ [A]®[B] ) = ( [B)®[A] )-[U 2m ] 

This matrix is frequently found in quantum information theory [Sj 
where one write, by using the Pauli matrices 



[U 2 m] = U o°®o° + <rW ) 



where, 
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are the Pauli matrices, a c 
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and according to Einstein convention 



repeated index denotes summation over the index in question. We call this 
matrix a tensor commutation matrix 2 x 2. A tensor commutation matrix 
3 x 3 has been written by KAZUYUKI FUJII |Sj by the following way 
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in order to obtain a conjecture of the form of a tensor commutation matrix 
n ® n, for any n G N*. 

We have studied the properties of such matrices in the more general case, 
according to the RAO EL IN A ANDRIAMBOLOLONA approach in linear 
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algebra and multilinear algebra In establishing firstly, the theorems on 
linear operators in intrinsic way, that is independently of the basis, and after 
that we demonstrate the analogous theorems for the matrices. 

Define [£7 n ® p ] as the tensor commutation matrix n<g>p, n, p e N*, whose 
elements are or 1. 

In this article we have given two manners to construct [f/ n ®p] f° r an Y n an d 
pGff*, and we have constructed a formula which allows us to construct the 
tensor permutation matrix [f/ni®n 2 ®...®n fc (c)], ( ni, n 2 , . . ., n k ) eW and a 
formula which gives us the expression of their elements. From ((TJ it is normal 
to think to what about the expression of [t^^] by using Gell-Mann matrices. 
But at first we are obliged to talk a bit about the definitions of the types 
of matrices, and after that we are going to expose the properties of tensor 
product. 

Define [I n ] as the n x n unit matrix. For the vectors and covectors we have 
used the RAOELINA ANDRIAMBOLOLONA's notations 0, with overlining 
for the vectors, x, and underlining for the covectors, ip. 

Throughout this article IK = R or C. 

1 Matrices. A generalization 

If the elements of a matrix are considered as the components of a second 
order tensor, we adopt the habitual notation for a matrix, without bracket 
inside. Whereas if the elements of the matrix are, for instance, considered as 
the components of sixth order tensor, three times covariant and three times 
contravariant, then we represent the matrix of the following way, for example 
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[M] = (.U™) 

hi 2 i3 = 111, 112, 121, 122, 211, 212, 221, 222, 311, 312, 321, 322 
row indices 

jihh = HI, H2, 121, 122, 211, 212, 221, 222 
column indices 

The first indices i\ and ji are the indices of the outside bracket which we 
call the first order bracket ; the second indices %2 and j2 are the indices of 
the next bracket which we call the second order brackets ; the third indices 
Z3 and js are the indices of the most interior bracket, of this example, which 
we call third order brackets. So, for instance, Mf^l = 5. 
If we delete the third order brackets, then the elements of the matrix [M] 
are considered as the components of a forth order tensor, twice contravariant 
and twice co variant. 

Consider a more general case [M] = wnere ^he elements of 

[M] are considered as the components of a 2/c-th order tensor, k times con- 
travariant and k times covariant. The first order bracket is the bracket of a 
ni x mi- dimensional matrix ; the second order brackets are the brackets of a 
n,2 x m2- dimensional matrices ; . . . ; the k-th brackets are the the brackets of 
the n k x m fc -dimensional matrices. [M] = ( 7t s ) 1 < s < nin2 ... nfc> x < t < mima ... rofc if 
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the elements of [M] are considered as the components a second order tensor, 
once contravariant and once covariant. Then, [H] 



s = n k n k _ x . . .ri3n2(zi-l)+n fc n fc _i . . .n 3 (i 2 -l) + - • ■+n k (i k - 1 -l) + i k (1.2) 



t = m fc m fc _i . . . m 3 m 2 (ji-l)+m k m k _ 1 . . . m 3 (j 2 -l)+. . ■+m k (j k _ 1 +j k (1.3) 



The elements of the matrix [N] = (N l k J ) 
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with interior brack- 



ets, can be considered as the components of a third order tensor, twice con- 
travariant and once covariant. Then, for example, (Aj 2 ) = 1. 



2 Tensor Product of Matrices 

Definition 2.1. Consider [A] = (A)) E M mxn (K), [B] = e M pxr (K). 
The matrix defined by 

" A}[B] ... Aj[B] ... Ai[B] 



[B] 



A\[B] ... Aj [B] 



a;[b] 



Af[B] ... Af[B] ... A™[B] 
obtained by suppressing the brackets after the multiplications by scalar, 
j4*[-B], is called the tensor product of the matrix [A] by the matrix [B]. 

[A]® [B] e M 

mpxnr (K) 

Theorem 2.1. Consider the linear operators A e C(£,T) , B G C(Q,7i). 

[A] is the matrix of A in the couple of basis ((ej)i <i < n , {fj)i <j <m), [B] the 
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one of B in ((gk)i <k <d (^z)i <i <p) ■ Then, [A]<S> [B] is the matrix of A®B in 
the couple of basis (B, B\), where 

B = (eT®^T, el®^2, . . . , el ® ^, ® ^T, ® ^2 , • . • ,ei®~g~ T , . . . , e^®g!, e^® 
Define 

B = (e7® gt) 1 < i < nl < fc < r = ((e^<j <„) <8> ((Klip < r ) 

^1 = (^®^)i^< mil <,<p= ((/i)l<i<m) ® ((>u)l<I<p) 

Let 

[A]®[ J B] = ( 7t 8 ) 1 < 8 < mi , jl < t < np 

where, 

according to the formulas (jl.2j) and (J1.3|) 

s =p(ii - 1) +i 2 
t = r(ji - 1) +32 

One can also define elements of tensor product of the matrix [A] by the 
matrix [B] by the the following way 

[A] ® [B] = (c££) = (4/^) 

(cf. for example [B] ) where, 
Zi?2 are row indices 
J1J2 are column indices. 

In the first notation the elements of [A]® [B] are considered as components 
of a second order tensor , whereas in the second one, the elements of [A]® 
[B] are considered as components of a forth order tensor. As an example, 
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[A] ® [B] = r 
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the result of the calculus for the second notation, whereas 
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[A] ® [B] 



2 



3 



4 







1 



1 2 



2 



-2 



3 



4 4 6 8 



the result of the calculus for the first notation. 



Proposition 2.2. • [A,]) ® ([5 2 ] • [A 2 ]) = ([B^ ® [S 2 ]) • ® [A 2 ]) 



/or any matrices [Bi], [Ai], [B2], [A2] if the habitual matricial products [Bi] ■ 
[Ai] and [B2] ■ [A2] are defined. 

Proposition 2.3. [I n ] ® [J m ] = [J nm ] 

Proposition 2.4. Tensor product of matrices is associative. 

Theorem 2.5. Consider ([Ai])i<i<nxm a basis ofM nxm (^), ([ B j])i<j< pxr a 
basis ofMp X r(K)- Then, ([Ai] ® [5 j ]) 1 <.< rixm l < j < pxr a basis ofM npxmr (K). 

3 Tensor permutation operators 

Definition 3.1. Consider k IK- vector spaces Si, S2, ■ ■ ., o~ a permutation 
on {1,2,..., fc}. Call a- tensor permutation operator a linear operator U a 
from Si® £ 2 ®. . .® S k to £ CT (i)® £ CT (2)®- • -® £<r(fc), £(£1® £ 2 ®. . .® £ fe , 
£<r(i)<8> £<r(2)<S>- • -® ^tr(fc) ), defined by 

£4(^1 ® ^2 <8> • • • ® Xjfc) = ^(T)" <8> ^(2)" <8> • • • <8> ^(fe)" 

for all Xi E Si, X2 E S 2 , . . ., x~k E S k . 

If n = 2, then say that U a is a tensor commutation operator. 

Remark 3.1. U a is an isomorphism and U^ 1 is a tensor permutation oper- 
ator too. 

Theorem 3.2. Consider k K- vector spaces Si, S 2 , ■ ■ S k , a a permutation 

on {1, 2, ... , k}, U a E C(Si® S 2 ® ■ ■ ■ ® S k , S a{ i) ® £ ct(2) ® . . . ® 

cr- tensor permutation operator. Then, for all <pi E C(Si), (j) 2 G C(S 2 ), ■ ■ ■, 



4>k E C(S k ) 

U a -((f)l ® 2 ® - . . ® (f) k ) = (<f> a (l) ® CT (2) ® • • • ® 4>a(k))-U a 
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Proof. 4>i <g 4> 2 ® ■ ■ ■ <S> 4>k G £(£i <8> £2 <S> ■ ■ • <8> £*)> thus 

^•(01 ® 02 ®- • -® 0fc)G £(£l ® £ 2 ® • • • ® £fc, £<7(1) ® £<x(2) ® ■ • • ® 
0ct(1) <S> 0(7(2) <g • • • <g 0(T(fe) G £(£ CT (1) ® £<r(2) ® . ■ ■ ® £<7(fc)), thus 
(0o-(l) ® 0<t(2) ®- • -® (t><T(k))-U a e C(£l <g £2 ®. ■ -® £fc, £ CT (1) ® £x(2) ®- • -® 
£<r(fc) ) 

If xl e £1, z£ e . . Xfc e £ fc , 

^■(01 <8> 02 <8>- • -<8> 0fc) (SI (g) (g . . . <g Xfc) = £M0l(zT) <g 02(^2) ®- • -® 0fc(^fc)] 
= 0<r(l) (^(T)") ® 0<x(2)(^(2)") ®- ■ -® 0<r(fe)(^(fc)) 

( Since t/^ is a tensor permutation operator) 

= (0<x(i) ® 0,7(2) o ■ ■ • ® 0<r(fc))(^(Ty ® ® . . . ® x^y) 

= (0,7(1) <g CT (2) 8> ■ ■ ■ <g 0(7(fc))-C4(^l ® ^2 ® • • • <8> Xjfc). □ 

Theorem 3.3. If U a is a a- tensor permutation operator, then its transpose 
Uj is a a' 1 -tensor permutation operator. 

Proof. Let us consider k K- vector spaces £±, £2, £k of finite dimen- 
sions, U a G C(£i® £ 2 <g . • • <g £k, £ CT (i) ® £ CT (2) <8> • • . <8> £(7(fe)) cr- tensor 
permutation operator. Then Uj G £( £ CT (i)* ® £ CT (2)* ® • • • ® £o-(k)*> £\ ® 
£2 ® ••• ® £fe*)- 

Let G £ CT(1) *, ^) G £,(2)*, . . G Xi G £1, G £2, G £ fe . 

J7 CT * - (^W g) ^(2) (g, . . . <g) £W) (xl <g %2~ ®...®x£) 

= g y,<K2) g _ _ _ ^ ^(fc) [E/^Zl <g> X2~ <g . . . ® Xi)\ 

( by definition of operator transpose [Hj) 

= y9 CT W g yW 2 ) g . . . g ^W (x^iy <g <g . . . <g x^y) 

= ^ (1) (^Ml))®^ (2) (^))^- • -®^ (fc) (^w) = ^ (1) (^0))^ (2) (^(2T) • • .^(x 
(because (x^y) are elements of K) 



= ^(xi) <g> <£ 2 (x 2 ) <g . • . <g f_ k (x k ) 

= (cp 1 <g> ^ 2 (g) . . . ® (f k ) (xl (g X2" (g ... (g x^) 



We have 

Uj ( (f aW (g) ip a{2) <g) . . . <g> y? CT(fc) ) (xl <g) Xa" ® . . . <g) J*) 

= <g> y? 2 (g) . . . <g> (p k ) (xT <g> x^ <g ... ® x£) 
for all Xi G x-j G £ 2 , • • •, xTT G 
Hence, 

Uj ■ ( y CT(1) ® <g> . . . <g> y CT(fc) ) = <y ® ^ 2 <g . . . ® □ 

4 Tensor permutation matrices 

Definition 4.1. Let us consider k finite dimensional IK- vector spaces Si, £ 2 , 
. . ., £k of dimensions ni, n 2 , . . ., rik respectively, U a G C(£\® 8 % (g . . . <S> 
£k, £o{i) ® £0(2) ® ■ ■ ■ ® £ a {k)) °~- tensor permutation operator. Let 
$1 = (eli ® 612 ® ■ ■ ■ ® eirn) a basis of £ 1 
B 2 = (ejjT <8> <g ... ® e^Tj) be a basis of £ 2 

B k = (ekl <%> ek~2 ® ••• ® ej^) be a basis of £ fc 

[J7 CT ] the matrix of [7 CT with respect to a couple of basis (£>i®£> 2 ® . . .<g> 
-Bo-(i)® $<r(2)® • • - ® &a(k) )• [C^r] is a square matrix of dimension n\ x n 2 x 
. . . x rife independent of the basis Si, £> 2 ,. . ., Call this matrix cr-tensor 
permutation matrix n\ ® n 2 ® . . . ® n^. 

Theorem 4.1. [U a ] a a -tensor permutation matrix n\®n 2 ® . . .®Uk if only 
if, for all [ax] G M niXl (K), [a 2 ] G M n2Xl (K),. . ., [a k ] G M nkXl (K) 

[U a ] ■ ([«i] ® [0-2] ® ■ ■ ■ ® [ak]) = [« CT (i)] ® K-(2)] ® • • • ® [a<r(fc)] 

Proof. " =>- " Let 577 G £1, cE^ G £ 2 , . . ., 57/7 G [«i], [a 2 ], ■ ■ ., unicolumn 
matrices formed by the components of 57T, 02,. . ., aTfc respectively with respect 
to B u B 2 ,..., B k . 

U a (a{ ®Wi <g ■ ■ ■ ®otk) = o~7jy 577^7 <g> . . . (g 577^ 

The components of oTT (g> 767^ (g> . . . (g 57/7 with respect to £>i®£> 2 ® ...(g) 
are the elements of the unicolumn matrix [ai] (g [a 2 ] (8> • • • <8> and the 
otalj) ® 57^7 ® . . . ® oT^y ones, with respect to <6 CT (i)<g> B ai2 )® . . .(g) 6^), 



9 



are the elements of the unicolumn matrix [a CT (i)] <8> [«cr(2)] <8> ■ ■ ■ <8> [aw*;)]- 
Hence, 

[U a ] ■ ® [a 2 ] ® ■ ■ ■ ® [a*]) = [a ff (i)] ® [a CT (2)] ® • • • ® [«<r(fc)] 
" <^= " Suppose that for all [an] e .M niX i(K), [a 2 ] £ .M n2X i(K),. . ., 
K] G X nfcXl (K) 

■ (M <8> [02] <8> • • • <8> [«fc]) = [a CT (i)] ® [a<r(2)] <8> • • • <8> [a<r(fc)] 
Let oi G £ 1, 02 G £2, • • ■, Ofe G and Si, # 2 ,. . ., B k be basis respectively of 
Si, £2, ■ ■ ., £fc where the components of Oi, 02,. . ., a^T form the unicolumn 
matrices [oci], [a 2 ],..., 

G £(£1 <g> £2 <S> • • • ® £fc, £ CT (i) <S> fo-(2) ® ■ • • <8> £a{k)) whose matrix with 
respect to (Bi®B 2 ® . . .® B k , B CT (i)® B^)® ■ ■ -® B a ( k ) ) is [U a \. Thus 

U a (W[ ® 02" ® • ■ • ® «fc) = «^(T)" <8> a^j <g> ■ ■ ■ <8> a^j 
This is true for all aT £ £1, £ £2, ...,«&£ £&. 

J7 CT is a cr- tensor permutation operator and [U a ] is a-tensor permutation 
matrix. □ 

Theorem 4.2. Lei [J7 CT ] be a -tensor 'permutation matrix n\ ® n 2 ® • • • ® n k . 
Then, for all square matrices [Ai], [A2],. . ., [Ah], of dimensions respectively 
ni, n 2 , . . ., n k 

[U a ] ■ ([Ai] ® [A 2 ] ® . . . ® [A fe ]J = ([A CT( i)] g) [A ct(2) ] <g> . . . <g> [A CT(fe) ]) • [f/ CT ] 

Proof. Let Ai £ £(£1), A 2 G £ (£ 2 ),..., A fc £ £(£&) whose matrices with 
respect to £>i, £> 2 ,. . ., S^ are respectively [Ai], [A 2 ],..., [A fc ]. Then [Ai] ® 
[A 2 ] (g) . . . ® [A fc ] is the matrix of Ai ® A 2 ® . . . ® A k with respect to Bi®B 2 ® 
. . .® B k . 

But, Ai®A 2 ®. . .®A k £ C(£i ® £ 2 ® . . . ®£ k ) and A^i^A^)®- ■ -®A a{k) £ 

C{£ a (l) ® £ CT ( 2 ) ® ••• ® £a(k)), thus 

^•(Ai®A 2 ®. . .®A fe ), (A^^A^)®. . .®A u{k) )-U a £ £(£1 ® £ 2 <g> . . . ® 

£k, £<t(1) ® £<t(2) ® ■■■ ® £a(k))- 

[Ar(i)] ® [At(2)] ® ■ ■ ■ ® [Ar(fc)] is the matrix of A (7( i ) ® A ct(2) ® . . . ® A a(k) with 
respect to B a{ i)® B a(2 )® ■ . .® B a(k) , thus ([A a( i)]®[A a{ 2)]®. . .®[A a(k) ])-[U a ] 
is the one of (A a m ® A a M ® ■ ■ ■ ® A a (k)) • U a with respect to (Bi®B 2 ® ■ ■ -® 
B k , B a{ i)® B a{2 )® ...® B a{k) ). 
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[U a ] ■ ([At] (g) [A 2 ] (g> . . . <g> [A k ]) is the matrix of JJ a ■ (A x ® A 2 <g> . . . ® A fc ) with 
respect to the same basis. 
Employing the Theore m^ .2^ 

U a ■ {Ay ®A 2 £ 

Hence 

[U a ] ■ ([AJ ® [A 2 ] ® 



. . . ® Afc) = (A CT(1) ® 4,(2) ® . . . ® • f/ CT 

. . ® [A fc ]) = ([At(1)] ® [4r(2)] ® ■ ■ ■ ® [At(*)]) • [^] 



□ 



Example 4.3. Let [a] 
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e M 3xl (K), \J3] 
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EM 2xl (K) 



[C4ig>2] ^ A^6x6(K) is a tensor commutation matrix 3 (g)2. Using the The- 
orei ri4.l( 

[*W( N®[/?] ) = [fl®H 
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(4.1) 
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is the tensor commutation matrix 2 ® 3. 

Using the Theorenl4~2l for all [A] G M 3x3 (K), [B] G M 2x2 (K) 
[U m2 }-([AMB]) = ([BMA})-[U m2 ] 

and 

[U 2 U([BMA}) = (lA}®[B})iU2® 3 } 

Remarks 4.4. 1. Tensor commutation matrices 1 <S> 12, 12 ® 1, 2 ® 6, 
6<g)2, 3®4et4®3 are 12-dimensional square matrices and are the 
only 12 x 12 square matrices which are tensor commutation matrices. 

2. For n G IN*, the tensor commutation matrix l£g>n is the n-dimensional 
unit matrix. 

3. For n G IN*, the tensor commutation matrix n (g> n is a symmetrical 
matrix. 

4. For n G IN*, for jo G IN*, tensor commutation matrices p ® n et 
n®p are inverses and transposes each other. 

5. If p is a prime integer number , the only p-dimensional tensor com- 
mutation matrix is the p-dimensional unit matrix, which is the tensor 
commutation matrix 1 ® p. 

6. Let E be a n-dimensional K-vector space, U G £(£ <8> £) a tensor 
commutation operator 

C7 = U' 1 

That is J7 is an involutive endomorphism. 

7. For n G IN*, for p G IN*, the set { [U mnp ], [U p ® n ], [U n ® p ] } equipped 
of the habitual product of matrices is a sous-group of the group of the 
permutation matrices np x np. 

8. We will remark { [Ui® 12 ], [U 2m ], [U m2 ], [U 3m ], [U 4 ® 3 ] } equipped of 
the habitual product of matrices is not a sous-group of the group of 
the permutation matrices 12 x 12 because [^/2®6]-[C4®3] is n °t a tensor 
commutation matrix. 

Notation 4.5. If U a is a tensor permutation matrix, then aU a ,(a G K) 
is a tensor permutation matrix too. So denote [^m®n 2 ®...®n fc (c)] the cr-tensor 
permutation matrix ri\ ® n 2 ® . . . ® formed by and 1. 
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Example 4.6. a is the permutation on {1, 2, 3}, such that cr(l) = 3, a (2) 
2, (7(3) = 1. 



For [a] ■ 
ren i4.1[ 



a 



of 



r 

r 



G A^2xi(K), using the Theo- 



[Uvgow^)] ■ (H ® IP] ® H) = [t]®[^]®N 

Using the same method as the one of the Examp 1 44.31 

1 



[^2®2®2(<7")] 



1 

1 

1 

1 

1 

1 

1 



5 Construction of a tensor permutation ma- 
trix 



For r e IN*, define 



as the matrix r x r whose elements are zeros 



except the z-th row and j-th column which is equal 1. The tensor commu- 



tation matrix n (g) p is a linear combination of some [E, 
1. 



ij J 



with coefficient 



E 



(n) 



<8) 



(p) 

A7 



(np) 

[p(i-l)+k}[ P (j-i)+l] 



(5.1; 



For example, the tensor commutation matrix [L^^] of (|4.1|) can be written 



[U 3 



E 



(6) 

li 



E. 



(6)" 
23 



E. 



(6) 
35 



+ 



E 



(6) 
42 



(6) 
54 



(6) 



or 

[^3® 2 ] 



E 



(2) 

11 



£0)1 
^11 



£7. 



(3) 
23 



(2) 
11 



(2) 
22 



(3) 
12 



+ 



E 



E- 



(3) 
32 



E. 



(2) 
22 



(2) 
12 



(3) 
21 



+ 



E. 



(2) 
21 



<8> 



-^33 



13 



We are going to determine the relation between i, j, k and I in order that 



the matrix 



E 



00 



E. 



(p) 



ki 



make up the linear combination giving the tensor 
commutation matrix n®p. Let [a] and [b] be two column matrices with n, p 
rows respectively. 





' a 1 ' 




' b 1 ' 


a] = 




Ab] = 






a n 




If 



Et } 



E 



(p) 

kl 



■ (H ® [b]) 



E 



(») 



a] I (8) 



E 



(p) 

kl 



lb] 



' 




' " 


a? 




b l 













a j b l 


where, in these three last matrices a\ b l , a?b l are respectively at i-th, 
/c-th, [p{i — 1) + A;]-th rows. 
On the other hand 

" b'a 1 



[b] ® [a] = 

where b l a? is at [n(l — 1) + j]-th row. 
Hence 



b l a j 



lfa n 
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\p(i-l) + k) = [n{l-l)+j] 



(5.2) 



(n) 



E 



(p) 



A.7 



the relation between i, j, k, I in order that E, 
terms of the sum giving the tensor commutation matrix n 
relation (|5.1j) . for each column 7 







p(j + l 



was among the 
p. Using the 

(5.3) 



the element at p{i — 1) + k — [nil — 1) + j]-th row is equal to 1 and the other 
elements of the same column are equals to zero. 



According to ()5.2|) and (|5.3|) the tensor commutation matrix n®p can be 
constructed by the following way without doing any calculus : 

Rule 5.1. Let us start in putting 1 at first row and first column, after that 
let us pass into second column in going down at the rate of n rows and put 
1 at this place, then pass into third column in going down at the rate of n 
rows and put l,and so on until there is only for us n — 1 rows for going down 
(then we have obtained as number of 1 : p). Then pass into the next column 
which is the (p + l)-th column, put 1 at the second row of this column and 
repeat the process until we have only n — 2 rows for going down (then we 
have obtained as number of 1 : 2p). After that pass into the next column 
which is the (2p + 2)-th column, put 1 at the third row of this column and 
repeat the process until we have only n — 3 rows for going down (then we 
have obtained as number of 1 : 3p). Continuing in this way we will have that 
the element at n x p-th row and n x p-th column is 1. The other elements 
are 0. 

Example 5.2. We have the tensor commutation matrix 3 <S> 5 
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[^3<8>5] — 



1 



























































1 



























































1 

















1 



























































1 



























































1 

















1 



























































1 



























































1 

















1 



























































1 



























































1 

















1 



























































1 



























































1 



Proposition 5.3. We have the following formula 



\U n ,i®n 2 ®...®n k { (T )\ 
(ni,...,n k ) 

E^{n 1 n 2 ...n k ) 
[n fc n. fe _i...n 3 n.2(ii-l)+n fc n fe _i...n 3 (i2-l)+...+^ 

(iiv>ifc)=( 1 >->i) 

Proof. For any [a r ] = K)^.^ G M nrXl , r = 1,2,...,*; 
[U ni ®n 2 ®...®n k (v)] ■ (M <8> M <8> • • • <E> [ojfc]) = [a CT (i)] (g) [a CT ( 2 )] ® . . . <g> [a a (k)] 
can be written as a sum of some matrices E^\ <S> E^- <S> ■ ■ ■ <S> E™ k - . 

*1.71 L 2J2 L kJk 

Determine the conditions on i t , j t , i 2 , J2, • • •, 4, jk so that E^ h ® E^ j2 <S> 
...(g) E^ jk make up the sum \U ni ®n*®...®n k (°)\- 

Jo-(l) Jcr(2) Jcr(fc) 



K(i)] ® K(2)] ® • • ■ ® K( fe) ] = (aj^aj^ . . . a a(k) 
with 



u l U 2 • • • a k — %-(!) U cr(2) • • • a a(k) 
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is found at [n^n^-i) ■ ■■ ?V(3)?V(2) (j CT (i)-l)+. 
1) + ^<r(fc)(i CT (fc-i) - 1) +i<7(fe)]-th row. 



•+?V(fc)?V(fc-l) • • • ™<t(3) ( j<r(2) — 



,.71 












4 2 


® . . . ® 










afaf 2 . . . a( fe 







where af, af 2 , . . ., a J k k and a-fa^ 2 . . . a J k k are respectively at ii-th, i2-th, . . -th, 
4-th and [n k n k -i . . .n 3 n2{i\-l)+n k n k - 1 . . .n 3 (i 2 -l) + . . .+n k (i k -i-i)+i k ]- 

, . ® E^- k is among the terms whose sum gives 



th rows. Thus El l h <g> 



[a)] if only if 



n k n k _ x . . . n 3 n 2 (ii - 1) + n^-i . . . n 3 (i 2 - 1) + . . . + n fe (z fe _i - 1) + i k 
= n^k)n^k-i) • • • n (T (3)n -(2)(j CT (i) - 1) + . . . + n -( fe )n -( fe _i) . . . ?V(3)Ut(2) - 1) 

+ ?V(fc)(icr(fc-l) - 1) + ja(fc)- 

However, 

~ [n k n k - 1 ...n3n2{ii-\)+n k n k -x...nz{i<2-l)+...+n k {^^ 

So for the column 7 

7 = n fc . . . n 3 n 2 (ji -l) + n k ... n 3 (j 2 - 1) + . . . + ra fc (j fc _i - 1) + j k 
the element at \n a {k)n a {k-x) . . . n CT ( 3 )n <T (2)(j -(i) -1)+- • •+n (T (fe)n (7 (fc_ 1 ) . . . n CT ( 3 ) (jo-( 2 ) - 
1) + n a{k)(ja(k-i) — 1) + jo-(fc)]-th row is equal to 1 and the other elements of 
the same column are zeros. □ 



17 



6 Expression of an element of a tensor per- 
mutation matrix 



Here n and p are any elements of N*. So it is a matter of generalizing 
the expression of an element of the tensor commutation matrix n ® n for any 
n G At first, study the above example for conjecturing the expression 

for the more general case. So we follow the way in the paper Then write 
[£^3<g>5] by the following way : 



3<g>5j 







i 

1 


(j 


(J 


u 


(J 




u 


u 


u 


u 


(J 




u 


u 


u 


(J 


(J 






n 


n 


n 


n 


n 




1 


n 


n 


n 

V.J 


n 




n 


n 


n 


() 


n 

V 7 








































1 


















" 


l 








" 




" 











" 




" 






































1 





































































1 















" 





1 





" 




" 











" 




" 









































1 





































































1 












" 








1 


" 




" 











" 




" 












































1 





































































1 









" 











1 " 




" 











" 




" 















































1 





































































1 


the rectangular matrices [I n xp] 






W) 


l<2<n,l<j<p 







(#0,^.^ ^ i where 5) is the Kronecker symbol. The matrix 

\ J I l<«<p,l<J<n J 

[/pxn] S [4x P ] = (S^) = 

where, 

iii 2 = 11, 12, . . . , In, 21,22, . . . ,2n, . . . ,pl,p2, . . . ,pn 
row indices, 

3x32 = 11, 12, . . . , lp, 21, 22, . . . , 2p, . . . , nl, n2,...,np 
column indices, 

is a n x p-dimensional square matrix, which suggest us the following propo- 
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sit ion. 

Proposition 6.1. 

Pmp] = - 

where, 

ilia = 11,12, In, 21,22,..., 2n, . 

row indices, 

Jlj2 = ll,12,...,lp,21,22,...,2p,. 
column indices. 



,pl,p2,...,pn 
, nl, n2, . . . ,np 



(6.1) 



Proof. Let [a] = (^)i^<„ e -Mnxi(K), [6] = (^ 2 ) 1 < j2 < p G AW (X) 



([«]®wr->(W-([<»]«i) 



«1«2 



= (M ® H) 



ft* 1 a 1 



1112 



□ 

Now, we are going to generalize the formula ()6.1|) . Let us consider the 



matrices 



<r(r) ' 



fe)l<i r <n CT(r) ,l<ir<nr' r 1,2, . . . , fc. 



[ n CT(1) xmj W ^n CT(2) xn 2 j W • • • W |^n CT(fc) xn fc j = (^X'.'jfe) = (^^2 " " " ^jl ) ' 

is a 72in 2 ... rife x n -(i)n CT (2) . . . n CT (fc) -dimensional square matrix, which suggest 
us the following proposition. 

Proposition 6.2. 

[U ni9na ®...® nk {<T)] = {UfX-X) = (^(1)^(2) ■■■ 5% L{k) 



Proof. For [a r 



a- 



Jr) 

r Jl<j r <n r 



E M nrXl (K), r = 1,2,..., k, 



([ai] <g> [a 2 ] ® . . . <g> [ajt]) 



2112. ..Ifc 



([^i (a)] • ([ai] (8) [a 2 ] . . . <g> [a fc ])) 



«l«2---«fe 



- rf 1 fi i2 rf fe n jl n j2 n jk 

" %i)%2) ' 2 " fe 

rii J<r(l)ri 2 „-M 2 ) X*fc „M*0 

i CT (l) a CT (l) J CT (2) a CT (2) • • • %fc) 

= a a(l) a a(2) ■ ■ ■ a a{k) 



(K(l)] ® K(2)] ® ■ ■ ■ ® Kw]) 
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□ 



7 Expression of [U^] by using the Gell-Mann 
matrices 

The Gell-Mann matrices [H] are 



Ai 



A« = -L 



i 



1 

1 


1 " 



1 
1 
1 
0-2 



A, 






-i 




1 
















i 





, A 3 = 







-1 







































-i ' 




" 





" 




" 














; A 6 — 








1 


, A 7 = 








— i 


i 










1 










% 






Ao denotes the 3-dimensional unit matrix. (Aj) 0<i<8 forms a basis of the 
vector spaces of the 3-dimensional square matrices. So according to the 
Theorem 12.51 the system (Ai <g> Aj) 0< i<8 o< j <s * s a basis of the set of 9- 
dimensional square matrices. Then, [^3^3] can be broken down into linear 
combination of this system. But being inspired by the expression of [^2^2] ( 
Cf. formula $IJ)) by using the Pauli matrices, we have directly tried to calcu- 
late the sum Ai ® A.;, i = 1, 2, . . . , 8. 



Ai ® Ai 



f A ® A 



100000000 
000100000 
000000100 
010000000 
000010000 
000000010 
001000000 
000001000 
000000001 



20 



Hence, 

[Unas] = |A ®A + \\ ® \ 

CONCLUSION 

We can construct a tensor commutation matrix, with or without calculus. 
We can also construct a tensor permutation matrix, but this time, by calculus, 
and we have the expression of an element of such matrix. So employing these 
matrices, a property of tensor product is in addition to we have already 
got. As the relation between tensor commutation matrix 2 ® 2 with Pauli 
matrices is frequently found in quantum Information theory, we hope that 
its analogous, the relation between tensor commutation matrix 3 ® 3 with 
the Gell-Mann matrices, will have also its applications in physics. 
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